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p ABSTRACT 

QJ • The influence of viscosity on the flow behaviour in sphericaUy symmetric accretion, 

I has been studied here. The governing equation chosen has been the Navier-Stokes 

■ equation. It has been found that at least for the transonic solution, viscosity acts as a 

CSJ I mechanism that detracts from the effectiveness of gravity. This has been conjectured 

to set up a limiting scale of length for gravity to bring about accretion, and the 
physical interpretation of such a length-scale has been compared with the conventional 
understanding of the so-called "accretion radius" for spherically symmetric accretion, 
iy-^ I For a perturbative presence of viscosity, it has also been pointed out that the critical 

Tij" . points for inflows and outflows are not identical, which is a consequence of the fact that 

^ sO [ under the Navier-Stokes prescription, there is a breakdown of the invariance of the 

stationary inflow and outflow solutions - an invariance that holds good under inviscid 
conditions. For inflows, the critical point gets shifted deeper within the gravitational 
potential well. Finally, a linear stability analysis of the stationary inflow solutions, 
under the influence of a perturbation that is in the nature of a standing wave, has 
indicated that the presence of viscosity induces greater stability in the system, than 
O ' has been seen for the case of inviscid spherically symmetric inflows. 
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1 INTRODUCTION 

The influence of viscosity, in a spherically symmetric accreting system, at least in qualitative terms, has been the object of 
the study that has been undertaken here. Accretion flows are governed by the equations of fluid dynamics. This in itself 
makes a convincing case for viscosity to be accorded some importance, since viscosity is an intrinsic property of fluids. Indeed, 
in accretion processes, viscosity does have a prominent role to play - as has been established particularly for the axially 
symmetric case of thin accretion discs, where viscosity, being purportedly a mechanism for the outward transport of angular 
momentum, effectively causes the infall of matter. 

However, the spheric ally symmetric a ccreting system has been studied very widely, using only inviscid hydrodynamical 
equations fBondi..l952: .Frank et al.lll992l] . In this study, instead of the inviscid Euler equation, the Navier-Stokes equation 
has been considered as one of the governing equations of the flow. F or a compressible viscous flow with constant coefficients 
of viscosity, the Navier-Stokes equation can be rendered in the form ibandau fc Lifshitzlll987i) . 

9v , VP GM, 1 I „2 /V A ^ ^'] 

_ + (vV) V + — + — r = - [^V\ + + C) V (V.v)J (1) 

For molecular viscosity, for which the actual magnitude may be quite small, the coefficie nts of viscosity rj and as 
has been mentioned before, have been taken to be positive constants Jbandau fc Lifshitzlll987i) . This greatly simplifies the 
mathematics, without largely compromising the underlying physics. In so far as viscosity is dependent on temperature (which 
can be connected to a dependence on the density of the accreting matter), its spatial variation is much less compared with 
that of the velocity of the flow, and in an actual accreting system where radiative processes facilitate efficient cooling and 
keep the sy stem somewhat close to being isothermal, it is an expectation that the assumption of constancy of 77 and largely 
holds true jBalbus fc Hawleylll998l) . 
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For the present purpose, 77 and have been tak en to be the conventiona l first and second coefficients of viscosity. They 
are expected to be of the same order of magnitude (iLandau fc Lifshit j|l987l) . Even as molecular viscosity is understood to 
be of not much significance in quantitative terms, its introduction in one of the governing equations should at least reveal in 
what qualitative manner would it influence the fiow dynamics - especially the behaviour of the critical points of the flows and 
the stability of the steady inflow solutions under the effect of a linearized perturbation. 

Considering now the vector identity, 

Vx (Vxv) = V(V-v) - V^v (2) 

it is seen that for an irrotational flow, such as in a spherically symmetric accreting system, the left-hand side of Q vanishes 
and leaves behind the condition, 

V (V-v) = (3) 

a result whose use could be made in Q. 

It is to be noted here that for a compressible fluid, the equation of continuity 

1^ + V- (pv) = (4) 

which is the other governing equation of the flow, does not give the result V-v = 0, as a steady-state solution. Comparing this 
fact with the steady solution of (0, and by using the identity given in @, it can be seen that the compressibility of accreting 
matter allows for viscosity to have a role to play, in a spherically symmetric flow under steady-state conditions. 



2 THE EQUATIONS OF THE FLOW 



For a spherically symmetric system, in which v = v(r,t) and p = p{r,t), the equations given by Q and are rendered 
respectively as. 



dt ^ dr p dr p \ 3 J dr 



J.2 Qj. 



and 

— -f — — ( vr'^') = 
dt dr ^ ' 

Use has also been made of a general polytropic equation of state 
P = kp"' 



(5) 

(6) 
(7) 



where 7 is the polytropic exponent with an admissible range given by 1 < 7 < 5/3, the restrictions being imposed by the 
isothermal limit and the adiabatic limit, respectively. 

The speed of sound Cs, with which the velocity of the flow is suitably scaled, is given by 



Cs = -ykp 



,7-1 



(8) 



The set of equations l|KJ-l|Hll, completely describes the problem. It can be seen that in the two coefficients of viscosity 
add up like simple scalars, which makes it possible to define a total viscosity 77tot = (4/3)?/ -f C,. 

In the steady state, explicit time dependence vanishes, i.e. (d/dt) =0, and from and @ will then follow the results 



dv 1 dP GM 
ar p dr r-' 



1 d 



J.2 flf. 



[yr^^ 



Vtot d 
p dr 

and 

respectively. From 1101 it further follows that 

1 d / 2\ V dp 

— -r \vr ) = T 

dr ^ ' p dr 

Using ||7||, JHJl and llll l. yields from @ the expression. 



d{^) 



+ /cff = 



2 dr 7 — 1 dr 
where the effective force /cff is given by. 



(9) 



(10) 



(11) 



(12) 
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fcB — 



GM 



rjtot d 



V dcs 
Cs dr 



(13) 



7 — 1 p dr 

The second expression on the right-hand side of 1131 is now to be studied closely. For accretion, it is a convention that 
V < 0. When r decreases, dcs/dr increases with a negative sign and for the transonic solution {—v/cs) increases with a positive 
sign. Hence, at least for the transonic solution, the product {—v/cs) {dCs/dr) increases with a negative sign for decreasing r. 



d_ 

dr 



f V dcs\ 
\ Cs dr ) 



> 



(14) 



On using this result in 1)13^ . it may be concluded that /off < fgj:!iv=GM/r^ . This implies that the presence of viscosity 
detracts from the effectiveness of gravity in spherically symmetric accretion. Since it is gravity that drives the process of 
accretion in spherical symmetry, viscosity has an inhibiting effect on it, which is quite contrary to what is seen in axially 
symmetri c accretion, where viscosity, in e ffecting the outward transport of angular momentum, actually aids the process of 



symmetri c accretion, wnere viscosity, m e n 
accretion JPringleHlQSlllFrank et alJll992l) . 



It can also be seen that with the introduction of viscosity in a governing hydrodynamical equation, it would be possible 
to obtain a limiting scale of length for the effectiveness of gravity to bring about the infall of matter. Viscosity is a microscopic 
phenomenon. Very far away from the accretor, the bulk influence of gravity, having been greatly weakened, becomes comparable 
with the inhibiting influence of viscosity at the microscopic level, and at such a distance, a scale of length may be obtained 
which, in a manner of speaking, may be understood to be a "viscous shielding radius", rvisc. By dimensional arguments it is 
also possible to have some quantitative understanding of TvIsc. The physical quantity rh/rjtot has the dimension of length and 
this quantity may be identified with ryjac- Here rh is the mass accretion rate, which is dependent on the mass of the accretor 
and the "ambient conditions" (jPrank et al.lll992l) . although for the viscous syst em being discussed here, this dependence - 
as opposed to that of the inviscid case, where the dependence is exactly known jFrank et alJll992^ - may be known only in 
an order-of-magnitude sense. With the identification that rviac = rn/rjtot, the "viscous shielding radius" would then conform 
to what can be intuited concerning the role of viscosity, that for a non-zero magnitude of r]tot, the length TvIsc would set 
a finite spatial limit to the accretion process. This understanding of the physical nature of rvisc can be compared with th e 
somewhat loosely defined concept of what has conventionally been known to be the "accretion radius" ijPrank et alJll993 l. 
given by rs,cc—GM/cs^ (oo). In all realistic cases with molecular viscosity, the condition rvisc S> ^acc prevails. The accreting 
system would then very likely be bound by the scale of length race. In such cases, the mass accretion rate rh, would be very 
close to the value that would be obtained for the inviscid case. 

Viscosity is also expected to infiuence the critical points of the flow solutions in the v'^ - r space. Critical points of the 
flow are given by the condition that in such a space, f or a first-order diffe rential equation in i; and r, both the numerator and 
the denominator of (i'^) /dr vanish simultaneously l|Chakrabartill99flll . The flow equation Q, however, is a second-order 
differential equation. So to get the critical point condition, ideally a first integral of Q is to be obtained, which, unfortunately, 
cannot be done for compressible flows. Hence an approximation is resorted to, that viscosity has a very small perturbative 
infiuence on the conditions governed by the inviscid equations. For molecular viscosity the validity of such an approximation 
is quite evident. Moreover, for large distances, the second derivative of v is comparatively small. The expectation here is that 
the critical point would be at a distance where this condition would hold. A cautionary note to be made here is that for the 
branch of the transonic solution very close to the accretor, the second derivative of v would be large. 

Using ||7|l, and llOI I. would then deliver from @, the result 



pr 



dv 
dr 



2cs 



GM 



rjtot 
P 



d^v 
dr^ 



2v 



(15) 



The relation p — p^o [cs/cs (oo)]^" is now to be used, where n = (7 — 1) ^ . It is then possible to obtain from the right 
hand side of (I15L the condition 



2cs 



GM 



+ 



rjtot 


Cs (00) 


2n 


(d\ 


2vc\ 


poo 


C-sc 




[^dr^ 


c rc2 ) 







(16) 



where the subscript label c, indicates critical point values. From 116|l . it is possible to have some idea of the influence of 
viscosity on the critical point of the transonic inflow solution. 
A rearrangement of terms in 1161 gives. 



GM ^ 7?tot 

2Csc^ poo 



{oc) 



1 _ i!^!^ 

2 Vc dr^ 



(17) 



For the transonic solution it may be said that a power law of the form v^r~ is followed. For free fall, 5 = 1/2. For the 
hydrodynamical case being discussed here, it may be supposed that in the vicinity of the critical point it will he 1/2 < S < 1. 
That will make the second term in parentheses in 11711 . a dimensionless positive number less than unity. Since for inflows, 
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convention has it that Vc < 0, it can then be said that the presence of viscosit y causes the critic al point to be shifted inwards 
along the r-axis, as compared with the inviscid case, where = GM/2csc' Iprank et alJll992l) . This, in itself, is consistent 
with the fact that viscosity has the effect of weakening the influence of gravity. At the critical point, the pressure term balances 
the gravity term. With gravity having been weakened, the pressure term is capable of balancing the gravity term deeper within 
the potential well. 

The use of the Navier-Stokes equation gives rise to another feature of the flow. A look at or its spherically symmetric 
adaptation @, makes it evident that unlike in the inviscid Euler equation, v makes an appearance here in the first power (in 
the viscous term). A direct consequence of this fact is that in the steady state, the Navier-Stokes equation governing the 
spherically symmetric flow, is not invariant under the transformation v — > — v. The occurrence of v in the first power in @, 
causes the invariance to break down. The important conclusion to be drawn here is that unlike in the case of inviscid flows, 
the critical points of the transonic inflow and outflow solutions, are n ot identical in th e - r space. It is straightforward 
to see that in the inviscid limit of 0, the invariance will be restored llChoudhurilll999l) and the conventional critical point 
conditions for the inviscid case will be obtained. 



3 LINEAR STABILITY ANALYSIS OF STATIONARY SOLUTIONS 

The steady state solutions of (|KJ and (|SJ are vp and po, on whi c h are imposed small flrst-order perturbations v' and p', 
respectively. Following a similar prescription by IPetterson et alj il980l) for inviscid spherically symmetric accretion, it is 
convenient to introduce a new variable, /=pur^, which in physical terms is the mass flux. Its steady-state value /o, as can be 
seen from © and itTITji . is a constant that is identified as the mass accretion rate. A linearized perturbation /', about /o is 
given by. 



/' = (^v'po + uop') 

Linearizing in the perturbation variables v' and p', gives from Q, the result. 



21 
dt 



]_df_ 

r'^ dr 



With the use of l|19|l . successive differentiation of IjlH^ with respect to time yields, first 



dv' 
dt 



I 



por- 
and then 



df 
dt 



- vo 



dr 



9V 

9*2 



1 

por2 



d 



_ 

df^ " dr 



df_ 

dt 



Now linearizing in terms of the perturbation variables, gives from the result 



dv' d , 2 p' 

-TTT + TT- VOV + CsO 

dt dr \ po 



Vtot 

po 



d_ 

dr 



]_d_ 

r^ dr 



(r^v'^ 



Po dr 



y.2 Qj. 



(r^vo) 



(18) 



(19) 



(20) 



(21) 



(22) 



where Cao is the speed of sound in the steady state. 

Partially differentiating I22II with respect to time, and using the conditions given by (EHl, ^ and ^ will finally deliver 
the perturbation equation as 



dt^ dr \ dt I Vo dr 



I 2 2\ df 

Vo [Vo - CsO -TT- 



2 I d 



r^ dr \ Po dt po dr 



+ ■ 



1 df d 



J.2 Qj. \ 



Vo) 



(23) 



por^ dr dr 

It is quite evident that the related perturbation equation for the inviscid case jTheuns fc Davidlll992ll . will be readily 
obtained for jjtot = 0. 

For 1231 1. a solution of the form /' = g{r)e^* is chosen, in which g{r) is the spatial part of the perturbation and Q in 
general is complex. This leads to a quadratic equation in Q, which is 



Bf7 C = 



(24) 



where 
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A 

B 



= g 



2±M 



rjtotr 



dr 



-2 dr 



dr 



Vtot < r — 



dr 



_l__d 

J.2 j^y. 



Vo_d£ 

Po dr 



1 dg d 
Po dr dr 



^±{r' 



(25) 



Stability of the subsonic inflows has been studied here with a particular emphasis. In such flows, it is readily understandable 
that for a perturbation which is in the form of a standing wave, there are two spatial points, one very close to the accretor 
and one very far awa,Y from it, where it is possible to have boundary conditions which would constrain the perturbation to die 
out ijPetterson et al.lll98(t ). Multiplying l|24|l by vog and integrating over the range of r, that is bounded by the two points 
where the perturbation dies out, will give. 



(26) 



where 
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B 



J Vog dr 

47r?7tot 
(-m) 

vo (no^ 
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■dr 
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Vtot 



Po 
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dr 



vq dpo 
po dr 



dr - 



- Vtot 



dr 



[yogr 



d 
dr 



Vo_d£ 

po dr 



dr (27) 



The results given by 12611 and 12711 have been arri ved at by introducing the simplification of partial integrations, and by 
requiring that all terms at the two boundaries vanish JPetterson et al.|ll98C[l. Use has also been made of an integrated form 



of IIUI I. which is povor 



(- 



i) /47r, for spherically symmetric inflows~ jFrank et alJll992^ . The integration constant 



physically identified with the mass infall rate and therefore m > 0. 
The solution of H26^ is given by. 



2A 



(28) 



B'^ _ C_ 
AA^ A 

The result in 1281 indicates that Q. has a real part and to ensure stability it is necessary that (fi) < 0. An examination 
of A and B will give. 



B_ 
A 



(^m) 



?;o g dr 



(vogr^) 



dr ; 



Anriu 



i-~vo)g dr 



[yogr 



dr 



(29) 



For inflows, the condition is (— wo) > 0, which, for a real g, would then imply that {B/A) > 0. Instabilities could still 
arise if {—C/ A) > 0, because Q, would then have one positive root. It is seen that 




In the inviscid term of 13011 . for subsonic inflows, the governing condition is uq^ < Cso^ ■ This would make the term always 
negative for any real g. For the viscous terms, the asymptotic properties of the integrands near the outer boundary, would 
have to be studied. For r — >cxd, the boundary conditions are po — 'poo, vo — >0 and g — >0. In that case, for the first viscous 
term, asymptotically at least, the products of the derivatives within the integrand, would be negative. For the second viscous 
term, the integrand in the numerator would be positive asymptotically, but having a negative denominator (since vo < 0), 
the term would be negative overall. It is to be expected that the conclusion drawn on the basis of the asymptotic argument, 
could be logically extended over the entire range of the integration. This would then lead to the conclusion that {—C/A) < 
for all situations of physical interest, and thus makes it possible to avoid having instabilities to develop in the system. 

It is t o be noted here that for inviscid subsonic inflows, a similar analysis would give the result that fl would be purely 
imaginary l|Petterson et aL. 198(1 '!. i.e. the perturbation would be oscillatory with a constant amplitude. The stability analysis 
for the viscous flow solutions developed above, indicates that even if the perturbation were to be oscillatory, its amplitude 
would decrease. Viscosity is fluid friction in its essence and as such it damps out the amplitude of the perturbation. This 
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would therefore establish that the presence of viscosity would induce greater stability in the system. This is a result that 
is somewhat surprising, in view of the fact that viscosity is a dissipative mechanism that releases heat into the system. To 
understand this seeming anomaly, it would be worthwhile to go back to ®, in which, the viscous term gives a mode for 
the loss of mechanical energy from the system. However, in the choice of equations here, no account has been taken of an 
equation for the energy balance, that, for viscous dissipation of the mechanical energy, would monitor the concomitant rise in 
the internal energy of the system. This loss of mechanical energy manifests itself as greater stability for the system, compared 
with the inviscid case. 

For the stability analysis of the transonic solution, a differ ent line of reasoning w ould have to be adopted, because there is 
no physically feasible inner bound ary condition for t he infl ow JPetterson et alJll980l) . and neither is the condition wo^ < Cso^ 
satisfled within the sonic radius. iTheuns fc DavidI (^O^) have shown that in the case of inviscid flows, for the subsonic 
branch of the transonic solution, the flow stability is depe ndent on cond itions prevailing at the outer boundary. For the 
supersonic branch, recourse is to be had to the argument of ICarlic^ lll979l) that for a finite supersonic region, a disturbance 
would be carried away in a finite time, and in this manner would ensure the stability of the solution. These arguments, it 
is to be expected, would hold good for the transonic solution for viscous fiows. The expectation is based on the analogous 
reasoning that for the inviscid fiow solutions, the conclusions drawn about the stability of the subsonic infiows are largely 
and qualitatively the same as those about the stability of the transonic infiow. More so, when, for viscous fiow solutions, it 
has been argued that the system shows itself to have greater stability for the subsonic stationary infiow solutions under the 
infiuence of a linearized perturbation. Hence, it may be safely established that a linear stability analysis indicates that in the 
absence of any energy balance equation, viscosity, as incorporated in the Navier-Stokes equation, would apparently induce 
greater stability in the stationary solutions in a spherically symmetric accreting system. 



4 CONCLUDING REMARKS 

It has been seen that the presence of viscosity in a spherically symmetric accreting system, adversely affects the effectiveness 
of gravitation in bringing about the infall of matter - at least for t he transonic s olution. To the extent that the transonic 
solution is the one that is most likely to be realizable jBondilllQsj : ICarlicklllQTgl : iRav fc Bhattachariejl20o3) . the effect of 
viscosity on such flows can have important implications - especially so, when a signiflcantly large and scale-dependent effective 
"turbulent viscosity" is taken into consideration. Such a prescription is not entirely without its context or relevance, because 
spherically symmetric infall of matter is expected to be seen in accretion from the interstellar medium on to an isolated 
accretor, and the interstellar medium, as h as been well known, displays turbulent behaviour. And since the Navier-Stokes 
equation is closely connected to turbulence l)Frischlll99fl^ . the extension to the use of an effective turbulent viscosity to study 
a spherically symmetric accreting system, acquires some justification for itself. 

Viscosity, as has been discussed in an earlier section, could possibly set up a limiting length-scale for gravity to be effective 
enough in driving a spherically symmetric accreting system. Yet it can be seen that for a weakly viscous system, the "viscous 
shielding radius" rviac(= rn/rjtot) would scarcely be able to set a reasonable spatial bound on the accretion process. In such 
a situation the accreting system would be spatially limited by the accretion radius race. However, the theoretical possibility 
is open that rvisc could become comparable with race with the introduction of a significantly large and scale-dependent 
turbulent viscosity, which could then make the viscous shielding radius define a noticeable limit on the active spatial range of 
the accretion process. The relevance of prescribing a scale-dependent turbulent viscosity becomes quite evident, when it is to 
be recalled once again that the interstellar medium, which lends itself readily to spherical accretion, is widely acknowledged 
to be a turbulent system. 

It has been seen here that if viscosity is considered to have a small perturbative infiuence on the conditions governed by 
inviscid equations, then the critical point for the infiow branch is marginally shifted deeper within the potential well. It is a 
matter worth some serious conjecture, as to how far reaching the consequences would be for the critical point, if instead of a 
small perturbative viscosity, a significantly large effective turbulent viscosity were to be imposed. Since, the transonic inflow 
solution passes through the critical point, its behaviour could get radically affected. 

Much the same thing can be surmised for the result of a linear stability analysis of the steady flow solutions, with a 
turbulent viscosity prescribed for the Navier-Stokes equation. The change could be even m ore significant, conside ring that 
such an effective turbulent viscosity, being scale-dependent, would have a spatial variation ijMeszaros fc Sill3ll977i ) - not to 
mention factoring in the accompanying rise in the internal energy of the system. 
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